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ABSTRACT 

We report on an X-ray observation of the Seyfert 1.5 galaxy Mrk 6 obtained with the EPIC instruments 
onboard XMM-Newton. Archival BeppoSAX PDS data from 18-120 keV were also used to constrain the 
underlying hard power-law continuum. The results from our spectral analyses generally favor a double 
partial-covering model, although other spectral models such as absorption by a mixture of partially 
ionized and neutral gas cannot be firmly ruled out. Our best-fitting model consists of a power law 
with a photon index of T = LBl^g' 2 , 2 , and partial covering with large column densities up to Ah ~ 
10 23 cm -2 . We also detect a narrow emission line consistent with Fe Ka fluorescence at 6.45]to'o4 keV 
with an equivalent width of 93^20 e V- Joint analyses of XMM-Newton, ASCA, and BeppoSAX data 
further provide evidence for both spectral variability (a factor of ~ 2 change in absorbing column) and 
absorption-corrected flux variations (by ~ 60%) during the ~ 4 year period probed by the observations. 

Subject headings: galaxies: individual (Mrk 6) — galaxies: Seyfert — X-rays: galaxies 



20) x 10 22 cm -2 ], likely due to the torus atmosphere. How- 
ever, detailed X-ray spectral modeling was limited by mod- 
est photon statistics and the fact that the absorption in 
Mrk 6 dominates the X-ray spectrum throughout most of 
the ASCA bandpass. 

We capitalized on the superior throughput (about an 
order of magnitude higher than ASCA) and good spec- 
tral resolution of the instruments onboard XMM-Newton 
to investigate further the X-ray emission and absorption 
properties of Mrk 6. We also used the XMM-Newton data 
to perform sensitive searches for temporal and spectral 
variability during the observation and in comparison with 
earlier ASCA and BeppoSAX observations. 

2. X-RAY OBSERVATIONS AND DATA REDUCTION 

2.1. XMM-Newton 

Mrk 6 was observed with the European Photon Imag- 
ing Camera (EPIC) onboard XMM-Newton (Jansen et al. 
2001) with on-source exposure times of 31.8 ks (EPIC p- 
n) and 30.8 ks (EPIC MOS 1 and MOS 2). It was also 
observed with the Reflection Grating Spectrometer (RGS) 
for 46.4 ks (RGS 1) and 37.9 ks (RGS 2). The EPIC 
data allow sensitive 0.2-12 keV imaging spectroscopy with 
moderate spectral (E/AE ~ 20-50) and angular (6" 
FWHM) resolution. The parameters of the individual ob- 
servations are listed in Tab. 1. 

Both the p-n and MOS data were acquired in 'full-frame 
mode' using the medium filter. The data were reduced us- 
ing the Science Analysis System (SAS) v5.3.3 with the 
latest calibration products. Due to background flares at 
the end of the observation period, 7% of the p-n and 31% 
of the MOS 1 and MOS 2 data were rejected. After inspec- 
tion of a source-free background spectrum extracted from 
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1. INTRODUCTION 

Markarian 6 (Mrk 6, IC 0450; z = 0.0188; 77 Mpc dis- 
tant for Hq = 75 km s~ 1 Mpc -1 ) is one of a handful of 
intermediate Seyferts of type 1.5, including NGC 4151 and 
Mrk 766, showing evidence for an 'ionization cone' rather 
than the 'ionization halo' expected for a type 1 viewing 
geometry (Meaburn, Whitehead, & Pedlar 1989; Kukula 
et al. 1996). A torus 'atmosphere' along the line of sight 
has been suggested to explain the presence of the ioniza- 
tion cone (e.g., Evans et al. 1993). The torus atmosphere 
is thought to be optically thick for ionizing radiation be- 
tween the Lyman edge and soft X-rays, but it allows ra- 
diation outside of this spectral range (e.g., broad optical 
line emission from the Broad Line Region) to pass through 
with only moderate obscuration. Mrk 6 shows significant 
optical line profile variations on time-scales of months to 
years, suggesting that at least some of the gaseous ma- 
terial creating the lines can undergo coherent variations 
(e.g., Rosenblatt et al. 1992; Eracleous & Halpern 1993). 

X-ray investigations of Seyfert galaxies are well suited to 
probing matter along the line of sight toward their central 
black hole regions. Prominent examples of intermediate 
Seyferts that have been well studied in X-rays include the 
'cousins' of Mrk 6, NGC 4151 (e.g., Weaver et al. 1994; 
Ogle et al. 2000; Yang, Wilson, & Ferruit 2001; Schurch & 
Warwick 2002) and Mrk 766 (e.g., Leighly et al. 1996; Page 
et al. 1999, 2001; Matt et al. 2000; Branduardi-Raymont 
et al. 2001). In contrast, the only published X-ray study 
of Mrk 6 is that of Feldmeier et al. (1999; hereafter F99) 
using 0.6-9.5 keV ASCA data. These data revealed heavy 
and complex X-ray absorption that was best fit by a double 
partial-covering model with large column densities [~ (3- 
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CCD 1 on the p-n and MOS detectors, we selected data 
only in the 0.3-12 keV energy range due to enhanced back- 
ground at low photon energies (up to a factor of <~ 5). We 
furthermore excluded bad and hot pixels from the data. 
Since the cleaned RGS data did not yield enough photon 
statistics for spectral analysis, they were not used further. 

Source counts were extracted for the p-n and MOS 
detectors within 30"-radius circular regions (the 85% 
encircled-energy radius at 7.5 keV). The background was 
extracted locally from source-free regions. For the p-n we 
used a rectangular region of size 2.'4 x 3f8 close to Mrk 6 
on the same p-n chip, and for the MOS we used an annu- 
lus with inner and outer radii of 2'5 and 4f3. Given the 
observed count rates of ~ 1.2 and ~ 0.4 cts s _1 for the p-n 
and MOS detectors, the pile-up fractions and dead times 
are estimated to be ^ 2%. The source counts were binned 
with a minimum of 50 counts per bin to allow x 2 spectral 
fitting. We also constructed background-corrected light 
curves of Mrk 6 in the soft (0.3-2 keV), hard (2-12 keV), 
and broad (0.3-12 keV) bands by binning the events into 
250 s intervals and excluding telemetry dropouts. 

2.2. BcppoSAX 

Archival BeppoSAX MECS 2+3 (Medium Energy Con- 
centrator Spectrometer) and LECS (Low Energy Con- 
centrator Spectrometer) data are also used in the anal- 
yses below (see Tab. 1); these data cover the 1.3-10 keV 
(MECS) and 0.5-4.5 keV (LECS) energy bands. To ex- 
tend the spectral range of our analyses to higher photon 
energies, which is useful to constrain the underlying power- 
law continuum, we furthermore used the BcppoSAX PDS 
(Phoswich Detection System) data in the energy range 18- 
120 keV. 

We checked for contamination of the PDS data due to 
other X-ray sources within its field of view (1?3 FWHM). 
No strong X-ray source other than Mrk 6 is found within 
the ~ 30'-diamctcr EPIC field of view in high-energy 
(> 5 keV) images. The ROSAT All-Sky Survey broad- 
band source catalog shows only two other X-ray sources 
(NGC 2256 and NGC 2258) located within the PDS field 
of view at large off-axis angles (23.'2 and 18^1) . Both 
of these sources have 2-10 keV MECS fluxes much lower 
(only ~ 2-3%) than that of Mrk 6. Since the PDS effective 
area decreases sharply with increasing off-axis angle, their 
contamination to the Mrk 6 PDS spectrum is negligible. 

3. RESULTS 

3.1. Spectral Analysis 

Spectral analysis was performed using XSPEC vll.2 
(Arnaud 1996). All errors are quoted at the 90% level of 
confidence for one parameter of interest (A% 2 = 2.71; Avni 
1976) and Galactic absorption (ATh = 6.4 x 10 20 cm~ 2 ; 
Stark ct al. 1992) is always implicitly included. The 
previous ASCA observation has shown that Mrk 6 is 
an intrinsically absorbed system (<; 10 22 -10 23 cm~ 2 ; 
F99). We therefore used the 18-120 keV BeppoSAX 
PDS data, which should be much less affected by ab- 
sorption, to determine the shape of the underlying X- 
ray continuum. Fitting a power-law model to the PDS 
data (model 1 in Tab. 2) gives a best-fit photon index 
of r = 1.81±°; 22 and an unabsorbed 20-100 keV flux of 
5.0 x 10~ n ergs cm~ 2 s _1 . 



We then fitted the XMM-Newton and BeppoSAX PDS 
spectra simultaneously. Due to uncertainties in the cross- 
calibration of the instruments (p-n and MOS: & 7%; Snow- 
den 2002) as well as possible variability between the XMM- 
Newton and BcppoSAX observations, we left all normal- 
izations for the spectral components free. All other spec- 
tral parameters were tied together. We first tried fitting a 
power-law model with simple (fully covering) intrinsic ab- 
sorption. Such absorption could arise in the host galaxy, 
perhaps in gas associated with the irregularly distributed 
dust seen in the HST image of Malkan, Gorjian, & Tarn 
(1998). It could also arise on smaller scales in the nuclear 
region. This model is statistically rejected (xt — 3.24; 

d. o.f. = 867). It also gives a photon index (r ~ 0.9) well 
below that derived from our BcppoSAX PDS analysis and 
expected intrinsically for a Seyfert galaxy (T ~ 1.6-2.2; 

e. g., Nandra et al. 1997; Risaliti 2002). 

We next tried fitting a power-law model with a single 
partial-covering absorption component. Partially covering 
absorption could arise from electron scattering of X-rays 
in the nucleus of Mrk 6 (see F99), and such absorption has 
been used successfully to explain the X-ray spectra of other 
intermediate Seyferts. This model also includes the sim- 
ple intrinsic absorption used in the previous model, since 
again there could be absorption in the host galaxy after 
the X-rays have escaped the nucleus (this simple intrinsic 
absorption will be implicitly included in all of the follow- 
ing spectral fits). Finally, we added a Gaussian emission 
line to model Fe Ka emission visible in the spectrum. This 
model is marginally acceptable (x 2 = 1.08; d.o.f. = 859) 
although it leaves clear systematic residuals below ~ 2 keV 
and above ~ 40 keV. Furthermore, like the previous model, 
it requires an implausibly low photon index of T ~ 1.3. If 
we constrain the photon index to lie in the range 1.6-2.2, 
the fit is statistically unacceptable. 

Given the failures of the two previous spectral models, 
we tried a model consisting of a power law and two par- 
tially covering absorbers (model 2 in Tab. 2). This model 
was successfully used by F99 to describe the ASCA data, 
and double partial covering could arise if electron scat- 
tering of X-rays provides several lines of sight into the 
nucleus. The photon index was fixed to the best-fit value 
derived from the BeppoSAX PDS data (r = 1.81; our 
results are not changed materially if the photon index is 
left as a free parameter), and we also included a Gaussian 
Fe Ka emission line. Despite some small systematic resid- 
uals below ~ 2 keV, a good overall fit to the XMM-Newton 
and BeppoSAX data is found (xl — 0.96; d.o.f. = 858; see 
Fig. 1). The absorption includes a high column density 

(10.96l°:42 x 1022 cm ~ 2 ) component covering (57 ± 1)% 
of the emitting region, as well as a lower column density 
(2.46t^o5 x 10 22 cm ~ 2 ) component with a higher covering 
percentage of (93 ± 1)%. The best-fit rest-frame energy of 
the Fe Ka line is 6.45±g;g| keV with an equivalent width 

of 93^20 eV- There is no clear evidence for a broad Fe Ka 
line, and the fitted line properties are consistent with those 
from ASCA. Given its properties, the observed line could 
originate via reprocessing in the outer part of the accre- 
tion disk (e.g., George & Fabian 1991) or perhaps in Broad 
Line Region clouds (see equation 5 in Eracleous, Halpern, 
& Livio 1996). 
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The Fe Ka line is fairly weak for a Seyfert galaxy, sug- 
gesting that Compton reflection is unlikely to make a dom- 
inant contribution to the X-ray spectrum. To examine any 
effects of Compton reflection further, we replaced the sim- 
ple power law in the previous fit with a power law plus 
neutral reflection (the 'PEXRAV model in XSPEC as- 
suming solar abundances). All other aspects of the model 
were unchanged. We obtain a relative reflection fraction 
(defined as r = Q/2ir, where fi is the angle subtended 
by the reflector) of r = 1.21^'^, an inclination angle 

of cos i = 0.94t°;°jt, and a cut-off energy of - 238 keV 
(model 3 in Tab. 2). The fit quality is slightly improved 
compared to the previous fit (xl = 0.95; d.o.f. = 855), but 
there is no material change in the nature of the partially 
covering absorption. Due to both statistical and system- 
atic uncertainties associated with this model, we do not 
consider the derived inclination angle to have clear physi- 
cal significance. 

We also tested a model which assumes that the nuclear 
X-ray emission is attenuated by both partially ionized and 
neutral material, dispersed along the line of sight. This 
model was proposed by Schurch & Warwick (2002) to char- 
acterize the complex X-ray absorption of NGC 4151, a 
Seyfert with many similarities to Mrk 6 (see §1). Although 
the model gives a relatively good fit to the 2-120 keV 
data, large residuals are visible below ~ 2 keV (x 2 = 1-14; 
d.o.f. = 866; see Fig. 1). We obtain an ionization param- 
eter of log £ = log (L ion /nr 2 ) = 2.42 ±0.01 (with L ion the 
source ionizing luminosity in the 0.0136-13.6 keV band, n 
the number of hydrogen atoms/ions per cm~ 3 , and r the 
distance from the central source to the photoionized re- 
gion), an ionized column density of 1.33+ ; 02 x 10 23 cm ~ 2 > 
and a neutral column density of 3.35io n x 10 21 cm~ 2 . 

A slightly better fit { X l = 1-09; d.o.f. = 863) is obtained 
if the simple power law in the previous fit is replaced with 
a power law plus neutral reflection (model 4 in Tab. 2). As 
for the double partial-covering model, the reflection is not 
a dominant spectral component and the best-fit absorption 
parameters do not change materially. 

We performed several additional fits to check the gen- 
eral robustness of our results. Models with no intrinsic 
absorption and Compton reflection cannot provide accept- 
able fits. Models with fully covering intrinsic absorption 
and Compton reflection can at best provide marginally ac- 
ceptable fits. These models also suffer from physical con- 
sistency problems since the energy region near the Fe Ka 
line is dominated by the reflection component. An Fe Ka 
line with a large equivalent width of ~ 1-2 keV is then 
expected (e.g., Matt, Brandt, & Fabian 1996) but not ob- 
served. 

Using the double partial-covering model (model 2 in 
Tab. 2), we derive observed 0.3-2 keV and 2-10 keV XMM- 
Newton fluxes of 5.3 x 10~ 13 ergs cm~ 2 s _1 and 1.2 x 
10~ n ergs cm~ 2 s _1 , respectively. These values are con- 
sistent with those in §2.3.5 of F99, given statistical and sys- 
tematic uncertainties. The absorption-corrected 2-10 keV 
flux and luminosity are 2.0 x 10~ n ergs cm~ 2 s _1 and 
1.4 x 10 43 ergs s _1 , respectively. The absorption-corrected 
0.3-2 keV flux and luminosity are less certain due to 
the large required absorption correction; likely values are 

5 The ASCA spectra used in this analysis were the same as those pi 



~ 1.5 x 10~ n ergs cm~ 2 s^ 1 and ~ 1.1 x 10 43 ergs s _1 , re- 
spectively. Using the partially ionized and neutral absorp- 
tion model (model 4 in Tab. 2), we derive an absorption- 
corrected 2-10 keV flux and luminosity ~ 36% lower than 
from the double partial-covering model. 

3.2. Temporal and Spectral Variability 

After fitting constant models to light curves and in- 
specting the residuals, we do not find any rapid, large- 
amplitude X-ray variability during the XMM-Newton or 
BeppoSAX observations. While constant-model fits are 
statistically inconsistent with the XMM-Newton p-n, Bep- 
poSAX MECS, and BeppoSAX LECS light curves, there 
are no strong, systematic variations about these fits (the 
poor fit quality appears to be due to stochastic small- 
amplitude variability only slightly larger than the statis- 
tical noise). We can constrain the amplitude of any sys- 
tematic variability to be 15% for the XMM-Newton p- 
n, <> 17% for the BeppoSAX MECS, and <, 29% for the 
BeppoSAX LECS. Within the individual observations, no 
evidence for spectral variability is found upon analysis of 
hardness ratios (e.g., 2-12 keV to 0.3-2 keV) computed as 
a function of time. No variability was seen in the ASCA 
data (F99), consistent with our results here. 

We have checked for inter-observation spectral varia- 
tions via joint spectral fitting. A joint spectral analysis 
of the XMM-Newton and ASCA data gives a fit with 
Xl = 1-04 (d.o.f. = 1595). 5 This joint fit is statisti- 
cally acceptable, indicating that we cannot prove spec- 
tral variability between the ASCA and XMM-Newton ob- 
servations. If we use the XMM-Newton and BeppoSAX 
data, however, we obtain xl = 1-20 (d.o.f. = 1410). This 
joint fit can be rejected with > 99.9% confidence, indi- 
cating spectral variability. Simultaneous XMM-Newton 
and BeppoSAX observations of the Seyfert 1 galaxies 
NGC 5548 (Pounds et al. 2003) and IC 4329A (Gon- 
doin et al. 2001) have shown that excellent fits can be 
found in joint spectral analyses, ruling out major cross- 
calibration errors. The spectral differences between the 
XMM-Newton and BeppoSAX data of Mrk 6 are largest 
below ~ 3 keV where absorption dominates the spectral 
shape (see Fig. 2), so absorption variability seems a likely 
explanation. We therefore analyzed the BeppoSAX data 
(MECS, LECS and PDS) individually, which gives best- 
fit column densities a factor of ~ 2 smaller compared 
to XMM-Newton (N^ = 4.761^ x 10 22 cm" 2 ; N^ = 

1.37lg'o7 x 10 22 cnl ~ 2 ) an d covering fractions which are 
consistent with our XMM-Newton results [/* = (54±1)%; 
/ 2 = (95 ± 1)%]. The observed 0.3-2 keV and 2-10 keV 
BeppoSAX fluxes are 1.1 x 10~ 12 ergs cm~ 2 s _1 and 2.4 x 
10~ n ergs cm~ 2 s _1 , respectively. These values are signif- 
icantly higher than during the XMM-Newton observation. 
The absorption-corrected 2-10 keV flux and luminosity 
are 3.2 x 10~ n ergs cm~ 2 s _1 and 2.1 x 10 43 ergs s _1 , re- 
spectively. The 2-10 keV luminosity of Mrk 6 during the 
BeppoSAX observation was ~ 60% higher than during the 
XMM-Newton observation. 

4. DISCUSSION AND SUMMARY 

The analyses above have substantially tightened the 
constraints upon the X-ray absorption and emission prop- 

3d in F99 (sec Tab. 1). 
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crtics of Mrk 6. The XMM-Newton spectrum has <~ 7 
times as many counts and a wider bandpass than the ear- 
lier ASCA spectrum, and the 18-120 keV BeppoSAX PDS 
data have provided the first reliable determination of the 
underlying X-ray continuum shape (a critical quantity for 
modeling of the X-ray absorption). The absorption mea- 
sured in the new data can be fit acceptably with the same 
basic double partial-covering model used to fit the ASCA 
data (F99), providing substantially improved support for 
the applicability of this model to Mrk 6. The measured 
column densities are large with iVn up to ~ 10 23 cm" 2 , so 
absorption controls the shape of the X-ray spectrum up to 
~ 6 keV. As discussed in §4.2.1 of F99, the absorption seen 
in X-rays is substantially larger than expected from the op- 
tical reddening; the X-ray absorbing material may be dust 
poor. The small 0.3-1 keV residuals left by the double 
partial-covering model may be due to <~ 5 x 10 40 ergs s _1 
emission from ionized gas in the nucleus or host galaxy. 
These residuals are consistent with X-ray emission from 
the Narrow Line Region, similar to that observed in NGC 
4151 (Ogle et al. 2000) and the Seyfert 2 galaxies NGC 
1068 (Ogle et al. 2003) and Mrk 3 (Sako et al. 2000). The 
residuals around 0.8-0.9 keV could be Fe L emission, but 
this spectral complexity is difficult to model due to the 
limited low-energy photon statistics of the EPIC instru- 
ments and the insufficient signal-to-noise ratio of the RGS 
spectra. 

Motivated by recent studies of NGC 4151 (Schurch & 
Warwick 2002), we also tried fitting the measured X-ray 
absorption with another physically plausible model con- 
sisting of both ionized and neutral columns of gas. This 
model provides an acceptable fit to the X-ray spectrum 
above 2 keV, but it leaves substantially larger residuals at 
lower energies than does the double partial-covering model 
(visible in the lower panel of Fig. 1). Again these residuals 
may be plausibly explained by emission from ionized gas 
in the nucleus or host galaxy. 

Three pieces of evidence suggest that, at least at high 
energies, our X-ray observations have penetrated all the 
way to the black hole region of Mrk 6: (1) the 18-120 keV 
continuum shape measured by the BeppoSAX PDS is con- 
sistent with those of Seyfert 1 galaxies (compare with §3.2 
of Malaguti et al. 1999), (2) the relatively small equivalent 
width of the Fe Ka line (EW = 93±2o eV ) indicates sig- 
nificant dilution of the reflection continuum at 6.4 keV by 
direct power-law emission from the black hole region, and 
(3) the relative 2-10 keV, [O in], and far-infrared luminosi- 
ties are consistent with those of Seyfert 1 galaxies (see §4.1 
of F99). Given these results, we can have confidence that 



the 2-10 keV (L 2 -io = 1-4 x 10 43 ergs s" 1 from XMM- 
Newton) and 20-100 keV (L20-100 = 3.5 x 10 43 ergs s" 1 
from BeppoSAX) luminosities we have found for Mrk 6 
represent those of the intrinsic X-ray continuum (rather 
than just scattered X-ray emission, for example). The 
X-ray luminosity of Mrk 6 is ~ 5 times the average X- 
ray luminosity of NGC 4151 (L 2 -io ~ 3 x 10 42 ergs s" 1 ; 
£20-100 ~lx 10 43 ergs s" 1 ; Yang et al. 2001; Schurch & 
Warwick 2002; Schurch 2002). Its lower X-ray flux arises 
only because it is ~ 5.8 times more distant. 

Our analyses have also revealed significant X-ray ab- 
sorption variability during the 1.5 yr between the Bep- 
poSAX and XMM-Newton observations. In the double 
partial-covering model, the column densities of both par- 
tial covering components drop by a factor of ~ 2 while the 
covering fractions do not change significantly. If the fitted 
partial covering in fact indicates that we have both direct, 
absorbed as well as electron-scattered, unabsorbed views 
of the nucleus, the observed changes seem physically plau- 
sible. The column densities along the direct line of sight 
could change due to gas motions in the torus atmosphere, 
while the scattered fraction would not change if scattering 
occurs on large spatial scales. The mechanism of absorp- 
tion variability must allow the column density to undergo 
large fractional changes; such changes could perhaps arise 
due to bulk rotation of a torus or Poisson fluctuations in 
a relatively small number of obscuring 'clouds' along the 
line of sight. Our discovery of X-ray absorption variability 
from Mrk 6 (as well as the observed variability amplitude 
and timescale) is generally consistent with the absorption 
variability seen from other absorbed Seyferts (e.g., Risaliti, 
Elvis, & Nicastro 2002; Schurch & Warwick 2002). 

Based on our results, further monitoring of the X-ray 
absorption variability of Mrk 6 is merited to refine under- 
standing of its amplitude and timescale. Such monitor- 
ing could be significantly complemented by coordinated 
optical spectroscopy, given the known optical line profile 
variability (see §1). High-quality ultraviolet spectroscopy 
of Mrk 6 would allow connections to be made between 
the X-ray absorption and any ultraviolet absorption. The 
archival IUE data on Mrk 6 suggest intrinsic Mg 11 ab- 
sorption, but higher quality data are required for a proper 
study. 

We thank L.A. Antonelli, M.C. Eracleous, and N.J. 
Schurch for helpful discussions. The project was supported 
by NASA grants NAG5-9940 and LTSA NAG5-8107 and 
NAG5-13035. 
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Table 1 

X-Ray Observations of Markarian 6 







Observation 




Exposure 


Observatory 


Instr. 


ID 


Date a 


(ks) 


ASCA 


CIS 


75041000 


1997-04-07 


42.1 


ASCA 


SIS 


75041000 


1997-04-07 


36.4 


BeppoSAX 


LECS 


51067001 


1999-09-14 


49.3 


BcppoSAX 


MECS 


51067001 


1999-09-14 


109.4 


BeppoSAX 


PDS 


51067001 


1999-09-14 


52.0 


XMM-Newton 


p-n 


0061540101 


2001-03-27 


31.8 


XMM-Newton 


MOS 1 


0061540101 


2001-03-27 


30.8 


XMM-Newton 


MOS 2 


0061540101 


2001-03-27 


30.8 


XMM-Newton 


RGS 1 


0061540101 


2001-03-27 


46.4 


XMM-Newton 


RGS 2 


0061540101 


2001-03-27 


37.9 



( a ) Start of observation. 



Table 2 

Spectral Fitting Parameters of Markarian 6 



Model Description and Data Used a 
(1) (2) (3) (4) 



Parameter PL PL+2PC+GA RE+2PC+GA RE+PI+GA 

(PDS) (PDS+p-n+MOS) (PDS+p-n+MOS) (PDS+p-n+MOS) 

Absorption, N H (10 22 cm- 2 ) 

Photon index, T 

Al (cts kcV" 1 cm" 2 s- 1 ) 

Partial covering, (10 22 cm -2 ) . 

Partial covering, (%) 

Partial covering, (10 22 cm -2 ) . 

Partial covering, / 2 (%) 

Photon index, T 

Cutoff energy (keV) 

Inclination angle, cos i 

Reflection fraction, r 

A2 (cts kcV- 1 cm" 2 s" 1 ) 

Warm absorption, N n (10 22 cm -2 ) 
Ionization parameter, log £ 

Line energy, E (keV) 

Line EW (eV) 

A3 (cts cm" 2 s- 1 ) 

X 2 /d.o.f. ( X i) 



Note. — All model fits include Galactic absorption by a column density of TVh = 6.4 X 10 20 cm 2 . 

( a )PL: power law; PC: partial covering; GA: Gaussian emission line; PI: photoionized absorption; RE: power law with Compton 
reflection by neutral material (PEXRAV). 



9.48 



i 01+0.22 
1 ' 81 -0.20 
+0.12 x 10 _2 



3.61/9 (0.40) 



0.11 



+0.01 
0.01 



0.11 



+0.01 
0.01 



1.81 (fixed) 



5.39 



+0.05 
-0.06 



x 10" 



10.96 



+0.58 



0.42 

57+ 1 



2.46 
93 



l 

+0.07 
-0.05 
+ 1 
1 



6.45 
93 



1.85 



+0.60 



+0.03 
0.04 
+26 
20 



8.05 
56 

2.27 
91 



+0.69 
-0.48 
+2 
-1 

+0.03 
-0.08 
+ 1 
1 



4.68 



1.81 (fixed) 

238t*f 

94+0 01 
u - a -0.25 
+0.14 
0.16 
+0.02 -, n - 
0.02 x ±u 



1.2V 



6.47 
87 



+0.03 
0.03 
+20 
26 



x 10" 



1 68+ ' 39 X 



io- 



824.2/858 (0.96) 812.8/855 (0.95) 



0.34 



+0.02 
0.01 



4.59 



1.81 (fixed) 
97Qc;+383 
i,OJ -2:)(i« 

95+ - 01 

1 OO+0.21 
J^-JJ-O.Ol 
+0.03 1Q - 

0.06 u 



13.32±°i 2 

2 42+ - 01 
z -^ z -o.oi 



6.45 
84 



1.67 



+ 0.41 



+0.04 
0.03 
+ 19 

26 



x 10~ 



937.3/863 (1.09) 
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Fig. 1 — XMM-Newton EPIC p-n (top), MOS (lower) and BcppoSAX PDS (right; > 18 keV) spectra of Mrk 6 in the 0.3-120 keV band. 
The EPIC counts were binned into energy channels with a signal-to-noise ratio of S/N > 10 for plotting purposes. The best- fit double 
partial-covering model is plotted as the histograms (see model 2 in Tab. 2). The fit residuals from this model are shown in the middle panel 
(in units of cr). The fit residuals from the partially ionized and neutral absorption model are given in the lower panel (see model 4 in Tab. 2). 
The positive residual at ~ 38 keV appears to be statistical; we do not find any evidence for a calibration error at this energy, and our results 
are not materially affected by the exclusion of this data point. The inset shows the (un-binncd) 5-8 keV EPIC p-n spectrum around the 
Fe Ka line. 




Energy (keV) 

Fig. 2.— XMM-Newton EPIC p-n (top; 0.3-12 keV), BeppoSAX MECS (middle; 1.3-10 keV) and LECS (lower; 0.5-4.5 keV) spectra 
of Mrk 6. Counts were binned into energy channels with a signal-to- noise ratio of S/N > 10 for plotting purposes. The best-fit double 
partial-covering model is plotted as the histograms (see model 2 in Tab. 2). The fit residuals are shown in the lower panel (in units of <r). 
Note the inability to fit the XMM-Newton and BcppoSAX data simultaneously with the double partial covering model. 



